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A new method of extracting diffractive parton distributions is presented which avoids the 
use of Regge theory ansatz and is in much closer relation with the factorisation theorem 
for diffractive hard processes. 



1 Introduction 



Diffractive parton distributions functions (DPDF's) [T] are essential ingredients in the under- 
standing and description of hard diffractive processes. The factorisation theorem for diffractive 
Deep Inelastic Scattering (DIS) [2j enables one to factorise the diffractive DIS cross-section 
into a long-distance contribution, parametrised by DPDF's, from a short-distance one, pertu- 
batively calculable. Although DPDF's encode non-perturbative effects of QCD dynamics and 
therefore must be extracted from data, their dependence on the factorisation scale is predicted 
by pQCD 1 . Moreover the short distance cross-section is the same as inclusive DIS [2J so that 
higher order corrections can be systematically accounted for. Due to the factorisation theo- 
rem, DPDF's are universal distributions in the context of diffractive DIS and diffractive dijet 
cross-sections are well described by next-to-leading order predictions based on DPDF's 3 . The 
commonly used approach [31 [H [5JH] to extract DPDF's is to assume proton vertex factorisation, 
i.e. that DPDF's can be factorised into a flux factor depending only on xp and t and a term 
depending only on the lepton variables (3 and Q 2 : 

f t D ((3,Q 2 ,x P ,t) = f P/P (x F ,t) if (/3,Q 2 ) + fn/p(x P ,t) jf (AQ 2 ) + ... 

Each term in the expansion, according to Regge theory, is supposed to give a dominant contri- 
bution in a given range of xp, the pomeron (IP) at low xp, the reggeon (JR) at higher value 
of Xp and so on. The flux factor fp/p (fp/p) can be interpreted as the probability that a 
pomeron (reggeon) with a given value of xp and t couples to the proton. This approach assumes 
an arbitrary truncation of the trajectory expansion and requires that parton distributions of 
each trajectory (ff, ff, ■ ■■) should be simultaneously extracted from data. It therefore in- 
troduces a large number of parameters in the fit and it is potentially biased by the choices of 
the flux factors. Although it has been proven to be supported by phenomenological analyses 
within HERA-I data precision, it is not routed in perturbative QCD and might be not entirely 
satisfactory with the expected precision increase of HERA-II data. 
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2 The new method 



The alternative method we propose is instead inspired by the factorisation theorem [2] for 
diffractive DIS itself. The latter states that factorisation holds at fixed values of xp and t so 
that the parton content described by ff is uniquely fixed by the kinematics of the outgoing 
proton and it is in principle different for different values of xp (and t, eventually). This idea is 
realised in practice by performing a series of pQCD fits at fixed values of xp with a common 
initial condition controlled by a set of parameters {pi}. This procedure guides us to infere the 
approximate dependence of parameters {pi} on xp allowing the construction of initial condition 
in the {/3, xp} space to be used in a global fit, without any further model dependent assumption. 
For the fits at fixed xp we choose the following singlet and gluon distributions at the arbitrary 
scale Qq\ 

PW,Qt) = A q p B « (1 - f3f« e-°^i , 
pg{P,Ql) = A g e~^i, 

which have four free parameters. We further assume that all lights quark distributions are 
equal to each other. The exponential dumping exponential factor allows more freedom in the 
variation of the parameters C q at large (3 and we choose the gluon distribution to be a simply 
a costant at Qq [J. Such distributions are then evolved with the QCDNUM17 [7 program within 
a fixed flavour number scheme to next-to-leading order accuracy. Heavy flavours contributions 
are taken into account in the general massive scheme. The convolution engine of QCDNUM17 is 
used to obtain and i^ 3 ' 3 ' structure functions at next-to-leading order which are then 

minimised against HI data [3]. In order to avoid the resonance region, a cut on the invariant 
mass of the hadronic system X is applied, M 2 > 4 GeV 2 . Fixed x p-Rt results are sensitive to 
the choice of the mininum Q 2 value of data to be included in the fits. The inclusion in the fits 
of data for which Q 2 < 8.5 GeV 2 in general worsens the \ 2 an d induce large fluctuation in the 
gluon distribution. This effect has been already noticed in Ref. [4] and avoided by including 
in the fit only data for which Q 2 > 8.5 GeV 2 . The same strategy will be adopted here. Good 
quality fits have been obtained with the common initial condition for all values of cc^-bins [5]. 
The dependence of the parameters (as returned by the fits at fixed xp) on xp is shown in 
Fig. pp. Red dots are the results from pQCD fits at fixed xp. The singlet normalisation A q 
behaves as an inverse power of xp. In order to improve the description at higher xp, however, 
an additional term is also included: 

A q {x P ) = A qfi {x P ) A ^ (1 - x P ) A ^ . 

The gluon normalisation is compatible with a single inverse power behaviour of the type: 

A g {xp) = A gfl {xp)^'- 1 . 

The coefficients B q and C q which control the /3-shape of the singlet distribution are well de- 
scribed by: 

B q (x P ) = B qfi + B q:1 xp, 

C q (Xp) = C q .o + C q .l Xp . 

The following generalised initial condition 
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Figure 1: Parameters as a function of xjp. Red dots are the results from pQCD fits at fixed xp. 
The grey line are best-fit prediction from Xjp-combined fit. The bands represent the propagation of 
experimental uncertainties by using the Hessian method [10] . 



PZ(/3,QIx p ) = A q (x P )f3 B ^ (l-(3f^^ e -^i, 
P g(P,Ql,x P ) = A g (xp) e~^i , 



is then used to perform a xjp-bin combined fit. The combined ht has nine free parameters. Fol- 
lowing the procedure described in Ref. [5] , to each systematic errors quoted in the experimental 
analysis is assigned a free systematic parameters which is then minimised in the fit along with 
theory parameters. As for the single- rrjp fits, only data points for which M 2 > 4 GeV 2 and 
Q 2 > 8.5 GeV 2 are included in the fit. The latter has an appreciable sensitivity on the scale Qq 
due to the relative stiffness of the initial condition. The choice of Qq is then optimised perform- 
ing a scan which gives the best x 2 value for Qq = 2.3 GeV 2 . The best fit returns ax 2 = 166 
for 182 degrees of freedom which is of comparable quality as the one presented in Ref. [J]. The 
initial condition allows the singlet and gluon normalisation, A q and A g respectively, to have a 
different power behaviour. It is therefore interesting to notice that if the condition A q ,i — A g .\ 
is enforced, this results in a global increase of the % 2 to 171 units for 183 degree of freedom. 
If one further neglects the a; jp-dependence of B q and C q by setting B q \ = C q .\ = the % 2 
increases to 188 units for 185 degree of freedom. This is an a posteriori confirmation that not 
only diffractive parton distributions change their magnitude versus xp but also that a modula- 
tion in their /3-shape (for the singlet, in this case) is necessary to better fit the data. The initial 
condition at Qq = 2.3 GeV 2 as a function of j3 for different values of xp are shown in Fig. [2]. 
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Figure 2: Singlet and gluon initial condition at Qq as a function of fi for different a;jp-values. The 
bands represent the propagation of experimental uncertainties by using the Hessian method [10] . 



3 Conclusions 

We have outlined a new method to extract diffractive PDF's inspired by the factorisation the- 
orem for diffractive DIS. From a series of pQCD fits at fixed xp we were able to infere the 
dependence of parameters on such a variable and this allowed us to construct a generalised ini- 
tial condition without assuming neither proton vertex factorisation nor the existence of a series 
of Regge trajectories. The best-fit returns a x 2 /d.o.f. close to unity, as the Regge-based pQCD 
fit of Ref. 0], but in our opinion the new procedure treats the non-perturbative a: ^-dependence 
of the cross-section in a controlled and less model dependent way and it might be capable (or 
even necessary) to fully exploit the expected improved precision of HERA-II data [111 [T2l [T3] . 
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